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We explore one-point and two-point statistics of the Navier-Stokes-a/3 regularization model at 
moderate Reynolds number (Re ~ 200) in homogeneous isotropic turbulence. The results are 
compared to the limit case of the Navier-Stokes-a regularization model and high resolution di- 
rect numerical simulation (DNS) results. After reviewing spectra of different energy norms of the 
Navier-Stokes-a and Navier-Stokes-a/3 model, we present probability density functions (PDFs) and 
normalized PDFs of the filtered and unfiltered velocity increments as well as longitudinal velocity 
structure functions of both regularization models and DNS results. We highlight differences in the 
statistical properties of the unfiltered and filtered velocity fields entering the governing equations 
of the regularization models and discuss the usability of both velocity fields for realistic flow pre- 
dictions. The filtered velocity field is found to have physically more viable PDFs and structure 
functions for the approximation of DNS results, whereas the unfiltered velocity field is found to 
have flatness factors close to DNS results. 

PACS numbers: 47.27.E-, 47.11.-j, 47.27.Gs, 47.27.ep 
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I. INTRODUCTION 



In principle, the accurate prediction of turbulent flows 
implies the resolution of the motion of the fluid on all 
spatial scales ranging from the large ones determined by 
the geometry of the flow domain to the smallest ones, at 
which the kinetic energy of fluid particles is transformed 
into heat by molecular dissipation. Even with access to 
state-of-the-art supercomputers, direct numerical simu- 
lations (DNS) of turbulent flows, meaning the solution 
of the Navier-Stokes (NS) equations at very high spatial 
and temporal resolution, are intractable, except for par- 
ticularly simple flow configurations. While DNS is widely 
used as a research tool to perform numerical experiments 
and explore the physical properties of turbulence, prac- 
tical flow predictions, e.g. in realistic engineering or geo- 
physical flows, are founded on turbulence models, which 
attempt to include the physical effects of turbulence at 
significantly lower computational costs than DNS. Tur- 
bulence models aim at predicting relevant large-scale flow 
features without resolving the smallest scales of the fluc- 
tuating turbulent velocity field. The interaction between 
the large scales and the small scales is parametrized in 
turbulence models. 

The class of regularized turbulence models stems orig- 
inally from the idea of performing statistical averaging 
at the level of the action principle using the Hamiltonian 
formalism and obtaining closure at the variational level 
with Taylor's "frozen-in turbulence" hypothesis, as dis- 
cussed by Holm The Navier-Stokes-a (NS-a) model 
(also called the viscous Camassa-Holm equations or La- 
grangian averaged Navier-Stokes (LANS-a) model) was 
first obtained by Chen et al. [2|-|4[ by addinga viscous 
term to the inviscid Camassa-Holm equations 5, 6], (also 
called the Lagrangian averaged Euler (LAE) equations or 
Euler-a equations). For a statistically homogeneous and 
isotropic flow with constant density p and constant dy- 



namic viscosity p, the NS-a model constitutes a system 

p(j^ + (gradv)u + (gradu) T v 

= -gradw + n(l - a 2 A)Au, } (LI) 

v = (l-a 2 A)u, divu = 0, 

for filtered and unfiltered velocities u and v related, as 
indicated in p.f j) 9. through the modified Helmholtz op- 
erator (I — a 2 A), and a filtered pressure-like variable 
vo = p-±p(|u| 2 +a 2 |D| 2 ), withD = ±(gradu+(gradu) T ) 
the filtered stretching tensor. In the limit a — » 0, the NS 
equations are recovered from 1|) . Various interpreta- 
tions of the filtered and unfiltered velocities have been 
provided by Holm et al. Q. In view of (jl.l)) ?. the ve- 
locity u is smoothed by applying the inverse modified 
Helmholtz operator (1 — a 2 A) -1 to the unfiltered ve- 
locity v, and thus contains less information than v. The 
constant filter parameter a has dimensions of length, and 
heuristically, represents the characteristic linear dimen- 
sion of the smallest eddy that the model can resolve. 
Numerical simulations indicate that the NS-a model cap- 
tures the properties of flows for eddy scales greater than 
a (see, for example, Chen et al. 



II. THE NAVIER— STOKES-q/? MODEL 

Based on a framework for fluid theories with higher- 
order gradient dependencies, Fried & Gurtin [jj-ITll] pro- 
posed a slight generalization of the NS-a model, which 
they call Navier-Stokes-a/3 (NS-a/3) model. For a sta- 
tistically homogeneous and isotropic turbulent flow of an 
incompressible fluid with constant mass density p and 
constant dynamic viscosity p, the NS-a/3 model consti- 
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tutcs a system 
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p(^jT + (gradv)u + (gradu) T v 



(1 



= -gradtn + fi(l - /J 2 A)Au, 
- a 2 A)u, divu = 0, 



(II.1) 

very similar to the system of the NS-a equations (jl.ll) . 
Whereas the parameter a > is the Helmholtz filtering 
radius, as in the NS-a model, and is meant to be repre- 
sentative of eddy scales in the inertial range, the param- 
eter /3 > is meant to be representative of eddy scales in 
the dissipation range; Fried & Gurtin [1-QjJ thus expect 
that j3 < a. Importantly, setting /3 = a in 1[) i yields 
the NS-a equations (|I.I[) and the NS equations are recov- 
ered in the limit a, /? — > 0. Notice that the only difference 
between the NS-a and NS-a/3 models is the modification 
of the viscous term as shown in (jll.ip i and (jl.ljh . 

In view of the identity grad(u • v) = (gradv)u + 
(gradu) T v + u x curlv and the relationship v = pip 
defining the kinematic viscosity, the flow equation (III. II) i 
can be written in the alternative form 

- u x curlv = -grad — + v{\ - /3 2 A)Au. (II.2) 
at p 



Multiplying (|II.2[) by u and integrating over the flow do- 
main R with periodic boundary conditions on dR yields 
an energy balance 



d 
di 



a 2 \u>\ 2 



)dv 



is(\u>\ 2 + /3 2 |grado;| 2 )d?;, (II.3) 



where u> = curlu is the filtered vorticity. With reference 
to (|II.3|) . the kinetic energy E and the dissipation rate e 
for the model, both measured per unit mass, are given 
by 



E = j i(|u| 2 +a 2 |u;| 2 )dw, 

R 

e = J v(\u\ 2 +^ 2 |grada;| 2 )dw. 

R 



(II.4) 



Notice that the kinetic energy (|II.4[) i of the NS-a/3 model 
is identical to that of the NS-a model. Setting /3 = a in 
(|II.4[) 9 recovers the dissipation rate of the NS-a model. 
Thus, choosing /3 < a reduces the dissipation of the ki- 
netic energy relative to the dissipation rate of the NS-a 
model. The equations in (|II.1[) are non-dimensionalized 
by introducing a characteristic length scale L, a charac- 
teristic velocity scale U, and a characteristic time scale 
L/U. For convenience, we maintain the symbols used for 
the original dimensional variables. As a result, the di- 
mcnsionlcss viscosity v represents the inverse Reynolds- 
number of the flow. The dimensionlcss length scales a 
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FIG. 1: Energy spectra E(k) of NS-a (a = 1/8) and NS-a/3 
(a = 1/8, /3 = 1/12) at 128 3 resolution, along with those of 
DNS at 256 3 resolution. 



and /3 arise simply by dividing their dimensional counter- 
parts by L. For numerical simulations of homogeneous, 
isotropic turbulence in a periodic box it is convenient to 
choose the dimensionless side length of the box to be 27r. 

Chen & Fried [l2[ used arguments similar to those of 
Foias et al. [HI to predict the scaling behavior of the NS- 
a/3 model in the inertial and dissipation ranges. They 
showed that, like the NS-a model, the NS-a/3 model is 
expected to recover Kolmogorov's fc -5 / 3 inertial range 
scaling and a faster k~ 3 falloff at higher wavenumbcrs. 
They offered a scaling analysis by which the wave number 
at which viscous dissipation becomes dominant is mod- 
ified by a factor a//3, implying that the inertial range 
of the NS-a/3 is extended by that factor with respect 
to the NS-a model. Kim et al. [U] studied the NS-a/3 
model numerically. They conducted numerical simula- 
tions of forced, homogeneous and isotropic turbulence 
in a periodic, cubic domain and showed that, provided 
reasonable choices of the two parameters a and /3, the 
NS-a/3 model, at lower resolutions, can closely approx- 
imate results obtained from highly-resolved simulations 
based on the NS equations. Confirming the theoretical 
predictions of Chen & Fried [l2| , they found that the en- 
ergy spectrum predicted by the NS-a/3 model for a > j3 
is more accurate than the energy spectrum from NS-a 
results at the same resolution. In recent numerical stud- 
ies Kim et al. [ll| investigated the utility of the NS-a/3 
model as a platform for spectral multigrid methods for 
periodic, homogeneous and isotropic turbulence. Sup- 
porting the view that the NS-a/3 model can approximate 
the NS-equations at lower resolutions, it was shown that 
the use of the NS-a/3 model at coarse grid levels and 
the NS equations at fine grid levels can accelerate the 
convergence rates compared to the convergence rates of 
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the same spectral multigrid method based on the NS- 
equations only. Kim et al. [HI investigated the usability 
of the NS-a/3 model for inhomogeneous flow with non- 
periodic boundary conditions. They presented a similar- 
ity theory for the NS-a/3 model to predict the microscale, 
which is found to be smaller than the microscale of the 
NS-a model, leading to the assumption that the NS-a/3 
model is able to capture smaller flow structures than the 
NS-a model, provided that a > /3. In numerical sim- 
ulations with a finite-element method, they considered 
the turbulent flow in two- and three-dimensional chan- 
nels past forward-backward steps. The results indicate 
that the physical features (in particular, eddy formation 
in the recirculation zone and eddy detachment) of this 
well-studied flow are qualitatively better captured by the 
NS-a/3 model for the choice of /3 = a/2 than the NS-a 
model at the same resolution and same a. These the- 
oretical and numerical studies of the NS-a/3 model sug- 
gest rough guidelines for choosing the length scales a and 
/3 < a both usually on the order of the mesh width. 



III. STATISTICS OF THE NS-a/3 MODEL 



Direct numerical simulations using (jTLTJ 2 .3 and (|II.2j) 
are performed for three-dimensional homogeneous and 
isotropic turbulent flows with periodic boundary condi- 
tions. We use the pseudospectral method for the spa- 
tial discretization and a second-order Adams-Bashforth 
scheme for time advancement. We employ a forcing 
scheme developed by Chen et al. [l7| that maintains con- 
stant energy in the first two wavenumber shells and, thus, 
extends the inertial range. The ratio between the energy 
contained in the first two wavenumber shells is chosen 
to be compatible with the fc -5 / 3 inertial range scaling. 
Unless otherwise indicated, all the simulations are per- 
formed using a (dimensionless) viscosity of v = 0.005 
corresponding to the Reynolds number Re = 200. 

All results are averaged over several large-eddy 
turnover times. We start collecting data once the flow 
reaches a statistically stationary state, which occurs af- 
ter approximately ten large-eddy turnover times. In 
the statistically stationary regime, we collect at least 
ten data sets over approximately thirteen large-eddy 
turnover times and then average those data sets. 

Numerical studies of regularization models available in 
the literature exhibit important differences on the strat- 
egy of comparison of results and evaluation of the quality 
of a certain model to produce approximations to DNS re- 
sults. We next briefly discuss the differences and describe 
our strategy. Most numerical studies of large-eddy sim- 
ulation (LES) models aim at comparing the predictions 
of a given LES model at a significantly lower resolution 
with a reference solution obtained from experiments or 
DNS data at a higher resolution. In most such stud- 
ies, results are not compared directly to DNS data at 
higher resolution, but instead a reference solution is usu- 
ally constructed by filtering high resolution DNS data 



and projecting the filtered data on the low resolution 
DNS grid, as explained by Meneveau & Katz [l8j]. For 
the evaluation and comparison of regularization models 
there has been examples of using unfiltered DNS data at 
higher resolution as reference solution. Chen et al. M, 
Mohseni et al. [l9| , and Pietarila Graham et al. [2(| EJ 
used this strategy in their numerical studies of the NS-a 
model. On the other side, Geurts et al. [22j adapted the 
LES strategy of using filtered and projected LES data in 
their numerical studies of different regularization models. 
The general LES strategy of resolving only the large, en- 
ergy containing scales while modeling the smaller scales 
including the dissipation range is based on the assump- 
tion that a separation of scales, normally only found in 
high or very high Reynolds number flows prevails. How- 
ever, work by Pietarila Graham et al.[2(|[2l| has demon- 
strated, that there are several difficulties with different 
regularization models at higher Reynolds numbers. Also, 
as opposed to LES, the NS-a/3 model aims to model a 
larger range of the energy spectrum, including part of 
the small scale motion in the dissipation range and, thus, 
comparing to filtered and projected DNS data at higher 
resolution is not suited to our objectives. In this study, 
we therefore adopt the strategy of comparing results of 
the NS-a and NS-a/3 model at a lower resolution to un- 
filtered DNS results at a higher spatial resolution. 



A. Energy spectra 

In this section, we study energy spectra for the NS-a/3 
model along with DNS data at higher resolution. Specif- 
ically, we highlight differences between the norms com- 
monly used for the discussion of results from regulariza- 
tion models, namely 
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E = - / v • udv, 
2 J 

R 

E u = lJ H 2 dv, E v = ^ / 

R R 



} (ni.i) 



Introducing the wavenumber k, we denote by E(k), 
E u (k), and E v (k) the respective spectra computed from 
the above norms. Since v = (I — a 2 A)u, (jlll.llh is 
equivalent to ([II.4[h and, thus, the only natural norm for 
the NS-a and NS-a/3 models, arising in the energy bal- 
ance (|II.3j) . Adopting LES terminology, the norm (jlll. 1 1) 
and the corresponding spectrum E u {k), involving only 
the filtered velocity u, can be interpreted as the spec- 
trum of the resolved kinetic energy. 

To begin, we display the natural energy spectra of the 
NS-a/3 and NS-a models at 128 3 resolution compared 
with DNS results at 256 3 resolution in Figure [T] Plots 
are provided for the NS-a (a = 1/8) and the NS-a/3 
(a = 1/8, /3 = 1/12) model. In the inertial range, the 
natural energy spectra for both NS-a and NS-a/3 follow 
Kolmogorov's —5/3 law and are close to the spectrum 
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obtained from DNS. Further, the natural energy spec- 
trum of the NS-a/3 model closely approximates the en- 
ergy spectrum of the DNS results in both the inertial 
and dissipation ranges. This result is consistent with the 
findings of Kim et al. [l4[. It indicates that, for j3 < a, 
the NS-a/3 model provides more energy at smaller scales 
than in the NS-a model. 

The alternative spectra of E u (k) and E v {k) of the NS- 
a and NS-a/3 models are shown in Figure [2] For both 
the NS-a and NS-a/3 models, the resolved kinetic energy 
spectrum E u (k) follows Kolmogorov's —5/3 law in the 
inertial range and is indistinguishable from the energy 
spectrum of the DNS results, as shown in panel (a) of 
Figure O However, the spectrum E u (k) of the NS-a and 
NS-a/3 models significantly deviates from the DNS en- 
ergy spectrum in the dissipation range. Throughout the 
higher wavenumbers, the resolved kinetic energy spec- 
trum of the NS-a/3 model is closer to the DNS energy 
spectrum than the corresponding resolved kinetic energy 
spectrum of the NS-a model. That is, even the filtered 
NS-a/3 velocity field contains more small-scale features 
than the filtered NS-a model. The energy spectra E v (k) 
obtained from the integral arc shown in panel 

(b) of Figure [2J Neither the NS-a nor the NS-a/3 results 
follow the fc~ 5 / 3 law in the inertial range. Further, both 
models significantly overprcdict the spectrum obtained 
from DNS. Strikingly, these inconsistencies are pervasive 
in the sense that they occur in both the inertial and the 
dissipation ranges. Although the NS-a results are closer 
to the DNS results than are the NS-a/3 results, neither of 
these models provides a good approximation to the DNS 
results. 

These results indicate that the natural energy spec- 
trum E(k) for the NS-a and NS-a/3 model shows good 
agreement with that of the DNS results. Importantly, 
the resolved kinetic energy spectrum of the filtered ve- 
locity fields provides a good approximation to the DNS 



results, with the NS-a/3 model showing better agreement 
than the NS-a model. In contrast, the spectrum of the 
unfiltered velocity field does not provide a suitable ap- 
proximation to the DNS results because of the significant 
overproduction of energy contained at intermediate and 
high wavenumbers. Consequently, whereas the filtered 
velocity field has physically more reasonable statistics in 
spectral space, the unfiltered velocity field exhibits un- 
physical statistics in spectral space. In the following, we 
confirm this observation with results from investigations 
of probability density functions (PDFs), structure func- 
tions, and flatness factors of the NS-a/3 and NS-a models. 

B. Two-point statistics 

Two-point correlations arc widely used to describe the 
spatial structure of turbulent velocity fields from a sta- 
tistical viewpoint. Here, we focus on a particular instant 
t of time and consider purely spatial correlations. For 
brevity and without loss of generality, we therefore sup- 
press dependence on t. Structure functions are two-point 
correlations, linking the velocities v(x) and v(x + r) at 
two different points x and x + r in space separated by a 
vector r through the velocity increment 

(5v(x,r) = v(x + r) - v(x). (III.2) 

The projection of the velocity increment <5v(x, r) onto 
the direction of r yields the scalar longitudinal velocity 
increment 

r 

5vi (x, r) = 5v(x, r) • f , r=- (III. 3) 

r 

where r = |r| is the separation distance with the di- 
mension of length. For statistically homogeneous and 
isotropic turbulence, the statistics of Svi are indepen- 
dent of the direction f but do depend on the modulus r. 
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FIG. 3: Same legend as Figure(2] Normalized PDFs P(Sui) and P(8vi) of the filtered velocity increment Sui and the unfiltered 
velocity increment 5v t for NS-a (a = 1/8) and NS-a/3 (a = 1/8, /3 = 1/12) at 128 3 resolution and DNS at 256 3 resolution with 
different normalized separation distances r — 5/32, r — 1/16, and r = 1/32. 



Hence, in homogeneous isotropic turbulence, the statis- 
tics of 5vi can be evaluated for arbitrary directions f . We 
non-dimensionalize the separation distance r with the do- 
main size L and keep the same symbols for convenience. 
The dimensionless separation distance r and the sepa- 
ration distance N r in mesh points are related through 
r = N r /N, with N being the number of grid points in a 
given direction. 

Let P(5vi) PDF of the unfiltered velocity increment 
Svi; similarly, let P((5vi) 3 ) denote the PDF of the third 
moment of the unfiltered velocity increment Svi- The 



PDFs P(Sui) and P((Sui) 3 ) of the filtered velocity incre- 
ments are defined analogously. To obtain better statis- 
tics, PDFs are evaluated with velocity increments in all 
three Cartesian coordinate directions x, y, and z and 
are averaged over ten data sets collected during approxi- 
mately thirteen large eddy turnover times. We compute 
PDFs of the NS-a/3 model for three values, r = 5/32, 
r = 1/16, and r = 1/32, of the normalized separation 
distances, corresponding, respectively, to separation dis- 
tances in mesh points of N r = 40, N r = 16, and N r = 8 
on the 256 3 grid. While large separation distances coin- 
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FIG. 4: Same legend as Figure(2] PDFs of the filtered velocity increment Sui and the unfiltered velocity increment Svi for NS-a 
(a = 1/8) and NS-a/3 (a = 1/8, /3 = 1/12) at 128 3 resolution and DNS at 256' ! resolution with different normalized separation 
distances r = 5/32, r — 1/16, and r = 1/32. 



cide with low wavenumbers in the incrtial range, small 
separation distances coincide with large wavenumbers in 
the dissipation range. 

Previous studies of DNS data by Chen et al. 



171 and 



Ishihara et al. [23} showed that the normalized PDFs at 
large separation distance are close to Gaussian. Further- 
more, the tails of the PDFs at smaller separation dis- 
tances are expected to exhibit nearly exponential behav- 
ior. This behavior is evident in the PDFs for our 256 3 
DNS results shown in Figure [3] along with the normalized 
PDFs of the velocity increments for the unfiltered and 



filtered velocity fields of the NS-cv and NS-cv/3 models at 
128 3 resolution. Consistent with the discussion of struc- 
ture functions provided in Section IIII C[ the PDFs are 
normalized with the square root of the second moment 
((Sui) 2 ) 1 ' 2 of the velocity increment at the corresponding 
separation distance r. 

The normalized PDFs of the NS-a and NS-a/3 models 
are indistinguishable from those of the DNS results in 
a neighbourhood of what appears to be a shared maxi- 
mum. However, the tails of the PDFs of the unfiltered 
and filtered velocity increment significantly deviate from 
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FIG. 5: Same legend as Figure [2] PDFs of the third-order moments of the filtered velocity increment Sui and the unfiltered 
velocity increment 5v t for NS-a (a = 1/8) and NS-a/3 (a = 1/8, j3 = 1/12) at 128 3 resolution and DNS at 256 3 resolution with 
different normalized separation distances r = 5/32, r — 1/16, and r = 1/32. 



those of the corresponding DNS results. In general, the 
normalized PDFs of the filtered velocity increments of 
the NS-a and the NS-a/3 model are narrower than the 
corresponding PDFs of the DNS velocity increments, as 
shown in panels (a)-(c) of Figure[3] On the contrary, the 
PDFs of the unfiltered velocity increments of the NS- 
a and NS-a/3 model have broader tails than the PDF 
from the DNS results, as shown in panels (d)-(f) of Fig- 
ure [3] This effect becomes more prominent with decreas- 
ing separation distance r. Since the normalized PDFs 
of the filtered and unfiltered velocity fields nearly coin- 



cide with the DNS results, the filtering operation has 
no influence on the large scales; however, the difference 
between the PDFs obtained from the filtered and unfil- 
tered velocity fields increases with decreasing separation 
distance r or, equivalently, with decreasing scale. Since 
small separation distances correspond to the dissipation 
range, we confirm the view that both regularization mod- 
els alter the small scale properties of the flow, while pro- 
viding a good approximation for the large scale proper- 
ties of the flow. For all three separation distances, the 
tails of the normalized PDFs of the NS-a/3 model are 
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slightly closer to the DNS results than those of the NS- 
a model. However, in both regularization models, the 
tails of the normalized PDFs of the unfiltered velocity 
field exhibit exponential behavior, even at intermediate 
separation distances. This result indicates that the unfil- 
tered velocity fields of both regularization models exhibit 
unphysical statistics, which improve to become more re- 
alistic through the filtering operation. 

Besides the discussion of the normalized PDFs in Fig- 
ure [H which gives an idea of the relative shape of the 
PDFs, the "raw" PDFs, shown in Figure 2] are useful to 
evaluate and compare the actual values of the filtered and 
unfiltered longitudinal velocity increments and to relate 
the real space statistics to the statistics in wavenumber 
space. Notice that the PDFs of the unfiltered longitudi- 
nal velocity increments of the NS-a and NS-a/3 models 
are, for all three choices of separation distance, signifi- 
cantly broader than the PDF of the velocity increments 
of the 256 3 DNS results. The NS-a/3 model gives higher 
probabilities of larger unfiltered velocity increments than 
the NS-a model at all separation distances, but the dif- 
ference between the PDFs of NS-a and NS-a/3 model 
increases with decreasing separation distances. This ob- 
servation seems consistent with the notion that /3 is as- 
sociated with the dissipation range. The PDFs of the 
unfiltered longitudinal velocity increments of the NS-a 
model are closer to the PDFs obtained from the 256 3 
DNS results than are the PDFs of the unfiltered longitu- 
dinal velocity increments of the NS-a/3 model. As already 
mentioned in the discussion of the normalized PDFs of 
the unfiltered longitudinal velocity increments, the un- 
filtered velocity fields of the NS-a and NS-a/3 model do 
not appear to possess physically good statistical proper- 
ties when compared to the statistical properties of DNS 
results at higher resolution. 

The PDFs of the filtered velocity increments are shown 
in panels (a)-(c) of Figure HI For each of the separation 
distances considered, the PDFs of the filtered velocity in- 
crements of the NS-a and NS-a/3 model are very similar 
and provide good approximations to the PDFs obtained 
from the DNS results. The PDFs of the filtered veloc- 
ity increments of the NS-a/3 model tend to be slightly 
closer to the PDFs obtained from the DNS results. At 
large separation distances, the PDFs of the filtered ve- 
locity increments of the NS-a and NS-a/3 model closely 
approximate the PDF of the DNS results and the differ- 
ence between PDFs of the NS-a and NS-a/3 model and 
the DNS results increases with decreasing separation dis- 
tance r. This confirms the view that large scale features 
of the flow are well approximated by the NS-a and NS-a/3 
model while smaller scales flow features are lost through 
filtering. An optimal choice of parameters a and /3 was 
not attempted in the current study and might well yield 
further improvements. 

Drawing a parallel between statistics in wavenumber 
space presented in the form of the energy spectra in Fig- 
ure [1] and [2] and statistics in physical space leads to the 
following explanation. The natural energy norm for the 



NS-a and NS-a/3 model involves the inner product v ■ u. 
The unfiltered velocity field v of the NS-a/3 model is 
observed to possess broader PDFs of the Svi velocity in- 
crements corresponding to higher probabilities of larger 
velocity increments, especially at intermdiate and small 
separation distances. Consequently, it appears logical 
that, at higher wavenumbers, the energy spectra of the 
NS-a/3 model are closer to the energy spectra of DNS 
than the energy spectra of the NS-a model. That is, 
through the availability of what might be viewed as an 
additional tuning parameter, namely /3 < a, the NS-a/3 
model allows for more small scale motion than the NS-a 
model. 

The PDFs P({Svi) 3 ) and P((Suif) of the third-order 
moments of the unfiltered and filtered longitudinal ve- 
locity increments are shown in Figure [5] They show 
the same qualitative behavior as the PDFs of the un- 
filtered and filtered longitudinal velocity increments, but 
the mentioned tendencies become more prominent in the 
third-order moments, representing a higher-order statis- 
tics of the velocity field. While the PDFs of the third- 
order moment of the unfiltered velocity field of the NS-a/3 
model and the NS-a model are very similar at large sepa- 
ration distances within the inertial range, the differences 
between the two models increase with decreasing separa- 
tion distance r, confirming the view that the parameter 
j3 is responsible for the properties of the model at small 
scales. From the PDFs of the third-order moments of the 
filtered longitudinal velocity increments at intermediate 
separation distance, we observe that the results from the 
NS-a/3 model are closer to the results obtained from DNS 
than the results from the NS-a model over a wide range 
of the PDF. 



C. Structure functions 

The scalar, longitudinal structure function of order n 
is defined as 

SM(r) = {Sv l (r) n ), (III A) 

where (■} denotes a suitable ensemble average. The 
"four-fifths" law as part of Kolmogorov's famous theory 
(see, for example Frisch [HI), dating from 1941 (and de- 
noted by K41 hereafter), relates the third-order velocity 
structure function in the inertial range of homogeneous 
isotropic turbulence with the dissipation rate e and pre- 
dicts a scaling 

S( 3 )(r) = {5vt(r) 3 ) = -\er (III.5) 

proportional to the separation distance r. More gener- 
ally, Kolmogorov concluded that the structure function 
of order n scales as 

SW(r) = {6vi(r) n ) = -C n (er)^ 3 , (111.6) 

where the coefficient C n is a universal constant. For 
higher-order structure functions, scaling exponents have 
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(d) (e) (f) 

FIG. 6: Second- and fourth-order structure functions of the NS-a (a = 1/8) and NS-a/3 (a = 1/8, f3 = 1/12) model at 128 3 
resolution compared to DNS results at 256 3 resolution and the K41 scaling prediction; (a)-(c) second-order structure functions 
of unfiltered, filtered, and mixed velocity increments (left to right); (d)-(f) fourth-order structure functions of unfiltered, filtered, 
and mixed velocity increments (left to right). 



been found to deviate from the K41 scaling prediction 
due to intermittency, as discussed by Frisch [24| and oth- 
ers. We discuss structure functions 

S<»>(r)HMr) n > (HI.7) 

of the unfiltered velocity field, structure functions 

Si n \r) = (5 Ul (rT) (111.8) 

of the filtered velocity field, and mixed structure func- 
tions 

Si n J(r) = ((Su l (r)6v l (r))^ 2 ). (111.9) 

In detail, the procedure we use to compute the structure 
functions in physical space is as follows. For a veloc- 
ity field at a fixed time i, we choose f to coincide with 
one of the Cartesian unit directions x, y, or z. For each 
separation parameter r, at each discrete grid point 
in the computational domain, we evaluate the velocity 
increments between the neighbouring points located at 
Xrf + ^x and x^ — ^x, respectively. The velocity incre- 
ments, or their higher-order moments are averaged over 
all points in the computational domain and determine 
one value (r) of the structure function. This proce- 
dure is repeated for all separation distances r in all three 
Cartesian directions. That is, we compute structure func- 
tions in the three Cartesian coordinate directions and 



average them to obtain more representative statistical 
results. This procedure is repeated for ten data sets col- 
lected during approximately thirteen large eddy turnover 
times and the results are averaged. 

Figure [5] shows the second- and fourth-order unfil- 
tered, filtered, and mixed structure functions of the NS-a 
(q = 1/8) and NS-a/3 (a = 1/8, (3 = 1/12) models at 
128 3 resolution compared to DNS at 256 3 resolution and 
the K41 scaling prediction. In the second- and fourth- 
order structure functions of the 256 3 DNS, we identify a 
short section around r w 0.15 where they scale with the 
corresponding K41 prediction. Notice that r w 0.15 cor- 
responds approximately to the same regime as r = 5/32, 
the separation distance for which we also computed PDFs 
of the longitudinal velocity increments presented and dis- 
cussed above. This confirms that r = 5/32 corresponds 
to a separation distance within the inertial range. The 
second- and fourth-order structure functions of the un- 
filtered velocity field of the NS-a and NS-cv/3 model ex- 
hibit significantly higher values than the 256 3 DNS along 
the full range of separation distances. The NS-a/3 model 
possesses higher moments of the unfiltered longitudinal 
velocity increments than the NS-a model. In the inertial 
range, the NS-a and NS-a/3 structure functions of the 
unfiltered velocity both appear to scale with a smaller 
exponent than the corresponding K41 scaling prediction. 
The second- and fourth-order filtered velocity structure 
functions of the NS-a and NS-a/3 model closely approxi- 
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mate the respective structure functions of the 256 3 DNS 
results. Near r « 0.15, the filtered structure functions 
of the NS-a and NS-a/3 model coincide in the inertial 
range; further, for smaller separation distances r, the 
NS-a/3 model exhibits higher values than than the NS- 
a model and, thus, in the dissipation range provides a 
better approximation to the 256 3 DNS results than the 
NS-a model. This behavior points to the presence of in- 
creased motion at smaller scales in the results of the NS- 
a/3 model compared to the NS-a model and is consistent 
with the behavior of the PDFs discussed in Section UlI Bl 
The mixed velocity structure functions of the NS-a and 
NS-a/3 model both show higher values than the corre- 
sponding 256 3 DNS. For separation distances r > 0.15 
the values of NS-a and NS-a/3 model coincide for both, 
second- and fourth-order mixed structure functions. At 
smaller separation distances, the NS-a/3 model displays 
higher values than the NS-a model and the mixed struc- 
ture functions of the NS-a model provide a better ap- 
proximation to the 256 3 DNS results. 



and NS-a/3 model. Interestingly, the unfiltered veloc- 
ity increments for both NS-a and NS-a/3 models exhibit 
slightly higher flatness factors than the velocity field from 
corresponding DNS results, as displayed in panel (a) of 
Figure [JJ This result indicates that the unfiltered ve- 
locity fields of those two regularization models are more 
intermittent than the DNS results at higher resolution. 
However, the flatness factor of the unfiltered velocity field 
of the NS-a/3 model closely approximates the flatness fac- 
tor of the DNS results over a wide range of intermediate 
separation distances, corresponding to the inertial range. 
That the unfiltered velocity field shows even higher flat- 
ness factors, and thus is more intermittent than DNS 
results, confirms our view that the statistical properties 
of the unfiltered velocity field are not sufficiently good to 
be considered a physically viable approximation to the 
velocity fields obtained from DNS results. 



IV. SUMMARY AND CONCLUSIONS 



D. Flatness of the velocity increments 

Central to the K41 theory is the assumption that the 
turbulent scales of motion are self-similar. This implies 
that the statistical features of turbulent flow should be 
independent of spatial scale. This self-similarity can be 
measured by the non-dimensional ratio of the fourth- 
order velocity structure function and the squared second- 
order velocity structure function. The flatnesses factors 
J- v and J- u of the unfiltered and filtered velocity incre- 
ments are defined as 



Mr) 



(Svfl 



Mr) 



(Suf) 



(111.10) 



respectively. If a flatness factor is flat or constant over 
all separation distances, the velocity distribution is self- 
similar. Otherwise, it is not self-similar, but instead it is 
characterized by intcrmittency. 

As observed in studies of DNS reported by Cao et 
al. [25jj, the flatness factor of NS turbulence is not con- 
stant as a function of the separation distance, due to the 
intermittent properties of the NS velocity fields. This 
observation is confirmed in our flatness results of our 
DNS results, shown in Figure [7] Chen et al. [1] reported 
flatness factors for the filtered velocity field of the NS-a 
model and found that they are systematically lower than 
those of corresponding DNS velocity fields. This indi- 
cates that the filtered velocity field of the NS-a model 
is less intermittent than the velocity field obtained from 
DNS results. In our simulation results, the flatness of the 
filtered velocity field of the NS-a/3 model is found to be 
very close to the flatness of the NS-a model, but slightly 
lower throughout all separation distances, as shown in 
panel (b) of Figure [7J Apart from the flatness factors of 
the filtered velocity field, here, we also discuss the flat- 
ness factors of the unfiltered velocity fields of both NS-a 



We presented and discussed Fourier and real space 
statistics of the NS-a/3 model in homogeneous isotropic 
turbulence. Findings were compared to the limiting case 
(/3 = a) of the NS-a model and to DNS results at higher 
resolution. We confirmed the result of Kim et al. fl4j 
showing that the natural energy spectrum of the NS-a/3 
can closely approximate DNS results at higher resolution, 
given a suitable choice of the parameter /3. We also found 
that the energy spectrum of the filtered velocity field, 
namely the resolved kinetic energy of the NS-a/3 model, 
contains more energy at the higher wavenumbers than 
the resolved kinetic energy of the NS-a model and, thus, 
provides a better approximation to high resolution DNS 
results. The energy spectra of the unfiltered velocity 
fields of NS-a and NS-a/3 model were both found to sig- 
nificantly overpredict the energy at intermediate and high 
wavenumbers. We found that choosing /3 consistent with 
/3 < a broadens the "raw" PDFs of the filtered and unfil- 
tered velocity increments, leading to higher probabilities 
of larger velocity increments. This corresponds to more 
small scale activity with ft < a and confirms the view 
that ft is responsible for the small scales. However, while 
the PDFs of the filtered velocity increments of the NS-a/3 
model slightly closer to the DNS results than the filtered 
velocity increments of the NS-a model, the PDFs of the 
unfiltered velocity increments of the NS-a/3 model were 
not found to provide good approximations to the DNS 
results. The structure functions of the unfiltered velocity 
field of the NS-a/3 model were seen to significantly over- 
predict the moments of the unfiltered velocity field, when 
compared to the structure functions of the DNS results. 
However, the structure functions of the filtered velocity 
field of the NS-a/3 model were found to exhibit better 
fidelity to the DNS results than the structure functions 
of the filtered velocity field of the NS-a model. For the 
mixed structure functions, the NS-a model showed bet- 
ter approximations to DNS results than did the NS-a/3 
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(a) 



(1)) 



FIG. 7: Flatness factors T v and T u of the unfiltered and filtered velocity structure function of the NS-a (a = 1/8) and NS-q/3 
(a = 1/8, P = 1/12) model at 128 3 resolution compared to 256 3 DNS results, (a) T v . (b) T u . 



model. As opposed to the structure functions, we found 
that the flatness, as a nondimcnsional, relative statistic 
of the unfiltered velocity field of the NS-a/3 model, shows 
excellent agreement with the unfiltered 256 3 DNS results 
over a wide range of separation distances, whereas the 
flatness of the filtered velocity field exhibits significantly 
lower values than the DNS results at small separation 
distances. These findings suggest that the unfiltered ve- 
locity field has intermittency properties very similar to 
the DNS results at higher resolutions, whereas the fil- 
tered velocity field is less intermittent. This is consistent 
with the idea of filtering. 

Our statistical studies indicate that choosing (3 such 
that f3 < a in the NS-a/3 model allows for more small 
scale activity in both the filtered and unfiltered velocity 
fields and, thus, provides good approximations to DNS 
results at higher resolution. Judging from the real space 



statistics of the NS-a/3 model, we conclude that the fil- 
tered velocity field has better, physically more reason- 
able statistical properties than docs the unfiltered veloc- 
ity field. The parameters a and j3 should therefore be 
adjusted in a way to aim for the best possible approx- 
imation of the filtered velocity field to DNS results at 
higher resolution. 
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